Unidirectional Invisibility induced by PT-Symmetric Periodic Structures 
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We show that parity-time (VT) symmetric Bragg periodic structures, near the spontaneous VT- 
symmetry breaking point, can act as unidirectional invisible media. In this regime, the reflection 
from one end is diminished while it is enhanced from the other. At the same time the transmission 
coefficient and phase, are indistinguishable from those expected in the absence of a grating. The 
phenomenon is robust even in the presence of Kerr non-linearities, and it can also effectively suppress 
optical bistabilities. 



In the last few years considerable research effort 
has been invested in developing artificial materials- 
appropriately engineered to display properties not found 
in nature. In the electromagnetic domain, such metama- 
terials make use of their structural composition, which in 
turn allows them to have complete access of all four quad- 
rants of the real e—fi plane. Several exotic effects ranging 
from negative refraction to superlensing and from nega- 
tive Doppler shift to reverse Cherenkov radiation can be 
envisioned in such systems [HE]- Quite recently, the pos- 
sibility of synthesizing a new family of artificial optical 
materials that instead rely on balanced gain/loss regions 
has been suggested [3HZj • This class of optical structures 
deliberately exploits notions of parity (P)and time (T) 
symmetry [8HT0] as a means to attain altogether new 
functionalities and optical characteristics [3]. Under VT 
symmetry, the creation and absorption of photons oc- 
curs in a controlled manner, so that the net loss or gain 
is zero. In optics, VT symmetry demands that the com- 
plex refractive index obeys the condition n(r) = n*(—r), 
in other words the real part of the refractive index should 
be an even function of position, whereas the imaginary 
part must be odd. PT-synthetic materials can exhibit 
several intriguing features. These include among others, 
power oscillations El HJ , absorption enhanced trans- 
mission [6J , double refraction and non-reciprocity of light 
propagation [3]. In the nonlinear domain, such pseudo- 
Hermitian nonreciprocal effects can be used to realize 
a new generation of on-chip isolators and circulators [7J. 
Other exciting results within the framework of PT-optics 
include the study of Bloch oscillations [12] , and the real- 
ization of coherent perfect laser absorbers [13] and non- 
linear switching structures |14j . 

To date, most of the studies on optical realizations 
of VT synthetic media have relied on the paraxial ap- 
proximation which maps the scalar wave equation to the 
Schrodinger equation, with the axial wavevector playing 
the role of energy. This formal analogy, allows one to in- 
vestigate experimentally fundamental VT -concepts that 
may impact several other areas, ranging from quantum 
field theory and mathematical physics [MIO], to solid 
state |T5] and atomic physics [TB]. Among the various 



themes that have fascinated researchers, is the existence 
of spontaneous VT symmetry breaking points (excep- 
tional points) where the eigenvalues of the effective non- 
Hermitian Hamiltonian describing the dynamics of these 
systems abruptly turn from real to complex 10J. Re- 
cently interest in VT -scattering configurations [T7H2U] 
has been revived in connection with using such devices 
under a dual role, that of a lasing and a perfect coherent 
absorbing cavity O . 




FIG. 1: (color online) Unidirectional invisibility of a VT- 
symmetric Bragg scatterer. The wave entering from the left 
(upper figure) does not recognize the existence of the periodic 
structure and goes through the sample entirely unaffected. 
On the other hand, a wave entering this same grating from 
the right (lower figure), experiences enhanced reflection. 

In this Letter we explore the possibility of synthesizing 
PT-symmetric objects which can become unidirection- 
ally invisible at the exceptional points. In recent years 
the subject of cloaking physics has attracted considerable 
interest, specifically in connection to transformation op- 
tics [TJ[55]. Here, our notion of invisibility stems from 
a fundamentally different process. As opposed to sur- 
rounding a scatterer with a cloak medium, in our case the 
invisibility arises because of spontaneous PT-symmetry 
breaking. This is accomplished via a judicious design 
that involves a combination of optical gain and loss re- 
gions and the process of index modulation. Specifically, 
we consider scattering from 'PT-synthetic Bragg struc- 
tures (see Fig. [I]) and investigate the consequences of VT 
symmetry in the scattering process. It is well known that 
passive gratings (involving no gain or loss) can act as high 
efficiency reflectors around the Bragg wavelength. In- 
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stead, we find that at the VT symmetric breaking point, 
the system is reflectionless over all frequencies around 
the Bragg resonance when light is incident from one side 
of the structure while from the other side its reflectivity 
is enhanced. Furthermore, we show that in this same 
regime the transmission phase vanishes-a necessary con- 
dition for evading detectability. Even more surprisingly, 
is the fact that these effects persist even in the presence 
of Kerr non-linearities. 

To demonstrate these effects we consider an optical 
periodic structure or grating having a PT-symmetric re- 
fractive index distribution n(z) = n + n\ cos(2/3z) + 
in 2 sin(2/3z) for \z\ < L/2. This grating is embedded 
in a homogeneous medium having a uniform refractive 
index hq for \z\ > L/2 (see Fig. [I]). Here m represents 
the peak real index contrast and n% the gain/loss periodic 
distribution. In practice, these amplitudes are small, e.g. 
n%,n2 C«o. The grating wavenumber /? is related to its 
spatial periodicity A via /3 = 7r/A and in the absence of 
any gain modulation (n 2 = 0) the periodic index modula- 
tion leads to a Bragg reflection close to the Bragg angular 
frequency = c/3/no (where c is the speed of light in 
vacuum). In this arrangement, a time- harmonic electric 
field of frequency uj obeys the Hclmholtz equation: 



d 2 E(z) u) 2 2 , , n 
— Y + -^n 2 (z)E(z) = 
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(1) 



For \z\ > L/2, Eq. M admits the solution E~(z) = 
EJ exp(zfcz) + EJ exp(-ikz) for z < —L/2 and Eq (z) — 
E~J exp(ifcz) + Ej cxp(-ikz) for z > L/2 where the 
wavevector k = hquj/c. The amplitudes of the forward 
and backward propagating waves outside of the grating 
domain are related through the transfer matrix M: 
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The transmission and reflection amplitudes for left (L) 
and right (R) incidence waves, can be obtained from the 
boundary conditions E£ = (EJ = 0) respectively, and 



are defined as t L = , r L = (t R = r R = ^+ , 

These can be expressed in terms of the transfer matrix 
elements as follows [T71 ITS] 
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While the transmission for left or right incidence is the 
same, this is not necessarily the case for the reflection. 
From the above relations one can deduce the form of the 
scattering matrix S [18] in terms of the M-matrix ele- 
ments. For "PT-symmetric systems, the eigenvalues of 
the iS-matrix either form pairs with reciprocal moduli or 
they are all unimodular. In the latter case the system is 
in the exact 'PT-phase while in the former one it is in the 



broken-symmetry phase [131 HZ] • For the complex peri- 
odic structure considered here, the transition from one 
phase to another (spontaneous PT-symmetry breaking 
point) takes place when n\ = n 2 |23j . 

To analyze this structure we decompose the electric 
field inside the scattering domain E(z), in terms of for- 
ward Ef(z) and backward Eb(z) traveling envelopes as 

E(z) = E f (z) exp(ikz) + E b (z) exp(—ikz). (4) 

Next by employing slowly varying envelopes for the field 
i.e. Ef(z) = E f(z) exp(iSz) and Eb(z) = £b(z) exp(— iSz), 
where S = j3 — k is the detuning. Substituting these ex- 
pressions in Eq. ([T]) , and keeping only synchronous terms 
while eliminating second order corrections in ni,and n 2 , 
we can then express the field at a point z inside the sam- 
ple in terms of the field at z = —L/2. For k f=s (3 close to 
the Bragg point, we get 
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where U = cos[A(z + L/2)]l - isin[A(z + L/2)]a ■ e, 
<7 are the Pauli matrices, and the unit vector e is de- 
fined as e = (1/A)(— fcn 2 /2no; — ikri\/2riQ\ S), while A = 
\J S 1 — k 2 (n i — n 2 , ) / Auq . By imposing continuity of the 
field at z = ±L/2, Eq. ^ becomes equivalent to Eq. ([2|. 
The transmission T = \t\ 2 and reflection coefficients 
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For n 2 = one recovers the standard scattering features 
of periodic Bragg structures. Namely, Rl = R R , while 
close to the Bragg point S — the reflection/transmission 
becomes unity/zero (in the large L-limit), see Fig. [2] In- 
stead if n 2 ^ 0, an "asymmetry" in the left /right reflec- 
tion coefficient starts to develop 23J. We would like to 
note that the VT arrangement considered here is funda- 
mentally different from that encountered in distributed 
feedback lasers (DFBs) [23]. In DFB systems both the 
index and gain/loss profile vary in phase and thus no 
■pT-symmetry breaking is possible. 

At ni = n 2 , this asymmetry becomes most pro- 
nounced. Even more surprising is the fact that at the 
Bragg point 6 = 0, the transmission is identically unity 
i.e. T = 1, while the reflection for left incident waves is 
Rl = (see Fig. |2|. This is a direct consequence of the 
VT- nature of this periodic structure. At the same time, 
the reflection for right incident waves grows with the size 
L of the sample as 
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of frequencies. For comparison, we also report in Fig. [3^,, 
(fit for the case of a passive (712 = 0) Bragg grating. 



8 (detuning) 



FIG. 2: (color online) Exact numerical evaluation (from Eq. 
|l])) of transmission T = \t\ 2 and reflection R — \r\' 2 co- 
efficients for a Bragg grating [25]. We have used no = 1, 
ni = 10~ 3 , L = I2.57T and /3 = 100. In case of a grating, 
the system is at the exceptional point when i%2 = n%. In this 



case, Rl is diminished (up to ni^ 



10" 



-see inset) for a 



broad frequency band, while Rr is enhanced. These exact 
numerical results are in excellent agreement with Eqs. ( 6|7 1. 



Such quadratic increase of the field intensity is a hall- 
mark of exceptional point dynamics [11) . This behavior 
is directly confirmed by our numerical simulations. We 
will refer to this phenomenon as unidirectional reflectiv- 
ity. Furthermore, Eqs. ^ indicate that a transformation 
ri2 — > — reverts the reflectivity of the system, allow- 
ing for reffectionless behavior for right incident waves i.e. 
Rji = 0, while the reflection from the left is now fol- 
lowing the prediction of Eq. tffh. In other words, the 
phase lag between real and imaginary refractive index 
dictates the unidirectional reflectivity of the system. 

Reffectionless potentials in one-dimensional scattering 
configurations are not in general invisible. This is due 
to the fact that the phase of the transmitted wave might 
depend on energy, thus leading to wavepacket distortion 
after the potential barrier. In this respect, a transparent 
potential can be detected from simple time-of-fiight mea- 
surements. It is therefore crucial to examine the phase (fit 
of the transmission amplitude t = \t\ exp(?0 t ) and com- 
pare it with the phase acquired by a wave propagating 
in a grating- free environment ((fit = 0)[53] Using Eq. ([5]), 
we deduce that the phase (fit close to the Bragg point is 
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At ni = ri2 we find that d — A, which results to a trans- 
mission phase (fit = 0. Thus interference measurements 
will fail to detect this periodic structure. Although the 
above theoretical analysis is performed close to the Bragg 
point 8 « 0, our numerical results reported in Fig. [3^,, 
indicate that these effects valid over a very broad range 




FIG. 3: (color online) (a) Transmission phase (fit as a func- 
tion of the detuning 8 for the PT-periodic system of Fig. [2] 
[25] . For comparison we plot together with the results of the 
"PT-exceptional point (ri2 = Hi), also the transmission phase 
for the passive structure (ji2 — 0). (b) The corresponding 
transmission delay times r t as a function of detuning. 

Next, we analyze the dependence of the transmission 
delay time Tt = d<fi t /dk [2H EZ], on the detuning 8. 
This quantity provides valuable information about the 
time delay (or advancement) experienced by a transmit- 
ted wavepacket when its average position is compared 
to the corresponding one in the absence of the scatter- 
ing medium. Using Eq. ^ we find that at the spon- 
taneous 'PT-symmetry breaking point the transmission 
delay time is r t = 0. In Fig. [3]d, we show results for a 
■pT-structure at n± = ri2 together with those expected 
from the passive case. 

It is also interesting to investigate the robustness of the 
above phenomena in the presence of Kerr nonlinearities. 
To this end, we assume the presence of a Kerr term in 
the refractive index profile i.e. n(z) = no + ni cos(2/3z) + 
in 2 sm(2f3z)+x\E(z)\ 2 . By decomposing the optical field 
into two counter-propagating waves and by considering 
only synchronous terms [23, [28] . we can then obtain a 
set of equations describing the field envelopes £b(z) and 
£f(z), in terms of Stokes variables [30] 

S (z) = 2kS 3 ; Si(z) = 2gS 3 ; S 2 (z) = 28S 3 - 3pS S 3 
S 3 (z) = -26S2 + 3pS a S 2 + 2KSo-2gS 1 (9) 

where p = k\/no, k = kn\/2no and g — knij2n§. 
It can be shown [23] that this non-linear system has 
the following conserved quantities gSo — kS\ = C%, 
3pgSQ—4:K8Si+4:KgS2 = C 2 - Using these constants of the 
motion, one can solve exactly Eqs. Because of lack 
of space we will not discuss the detail derivations here 
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[53] but rather cite the final results for the transmission 
and reflection coefficients 

(« + 0)So(f)-Ci D (k- 9 )S (-^)+C 1 

■ Hr — 



from AFOSR No. FA 9550-10-1-0433 and by the US- 
Israel Binational Science Foundation, Jerusalem, Israel. 
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In contrast to the linear case, now T L ^= Tr for n\ n<i 
indicating a diode action [7J [55] (see Fig. [4]) However, of 
interest here is the behavior of the system at the excep- 
tional point ri\ = ri2- We find that T L = Tr = 1, while 
Rl = 0, as in the linear case. These results are valid for 
any input intensity as shown in Fig. [4] At the same time 
we have found that the transmission phase is again inde- 
pendent of the detuning 8 and equal to <pt = 0. We thus 
conclude that the phenomenon of unidirectional invisibil- 
ity of "PT-periodic system at the exceptional point is en- 
tirely unaffected by the presence of Kerr non-linearities. 
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FIG. 4: (color online) Transmission and reflection coefficients 
in the nonlinear regime vs. left (left column) and right (right 
column) input intensities for 8 — (the same behaviour is 
observed for other values of 8). The parameters used are 
n = 1, ni = 0.5, and L = 7. Three different values of 
ri2 (below, above and at the exceptional point) are used. If 
712 7^ ni one can observe the standard bistability behavior of 
non -linear media. For ni = m, the bistability disappears, 
signifying the appearance of the spontaneous VT symmetric 
breaking point. At this point Rl = 0. 

We have shown that the interplay of Bragg scattering 
and VT symmetry allows for unidirectional invisibility 
which can be observed over a broad range of frequencies 
around the Bragg point. This process was found to be 
robust against perturbations. In the presence of nonlin- 
earities this unidirectional invisibility still persists and 
non- reciprocal transmission is possible. Of interest will 
be to investigate if these phenomena can also occur in 
higher dimensions and under vectorial conditions. 
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The Stokes variables are defined as So = + 
l^l 2 ,^! = \£ f \ 2 -\£ b \ 2 ,S 2 =£ f £i +£}£ b , S s = i(£}£ b - 
£ } £ b ) and they satisfy the identity S\ + S% + S$ = S 2 . 



